Evaluation of the economic characteristics of the fruit of 45 superior Camellia weiningensis Y.K. Li. trees

Reports related to Camellia weiningensis Y.K. Li. are rare. We evaluated the economic characteristics of the mature fruit of 45 superior C. weiningensis trees using principal component analysis (PCA) and gray correlation analysis, and identified excellent germplasms according to performance. PCA was employed to reduce the dimensions. PCA was performed for the original 15 indices of fruit diameter, fruit length, fruit shape, single-fruit weight, pericarp thickness, oil yield, fresh seed rate, dry seed rate, dry kernel rate and palmitic acid, stearic acid, linolenic acid, oleic acid, linoleic acid and arachidonic acid contents. According to the requirements of a cumulative contribution rate ≥ 80% and an eigenvector value > 1, six principle components were selected. These indices underwent weighted summation to establish a function model for comprehensive evaluation. Finally, the comprehensive rankings of the cultivars according to PCA were compared with those according to gray correlation analysis. The genetic variation coefficients of the 15 parameters ranged from 2.24% (oleic acid content) to 22.70% (single-fruit weight, with a range of 21.34 g). The top ten excellent cultivars with the highest comprehensive scores according to PCA and those according to gray correlation analysis were compared. According to PCA, oleic acid content, fruit diameter, fruit length, pericarp thickness, arachidonic acid content and dry seed rate can serve as representative evaluation indicators of C. weiningensis. The outcomes obtained based on PCA were basically consistent with those obtained based on gray correlation analysis. Finally, nine excellent cultivars were finally determined, i.e., WY-1, WY-6, WY-8, WY-25, WY-27, WY-30, WY-33, WY-35, WY-38 and WY-44. The results obtained in terms of crown yield were basically consistent with the outcomes of the comprehensive assessments, which indicates the reliability of the assessment methods used in this study.


Introduction
Camellia weiningensis Y.K. Li., belonging to the Camellia L. genus of the Theaceae family, is a type of evergreen shrub or small tree. The species is endemic to regions of 1800-2700 m altitude in Weining County, Guizhou superior cultivars have not yet been selected. In recent years, our research group has conducted systematic research on C. weiningensis resources and has collected and stored more than 200 excellent germplasms. After years of observations and comparisons, 45 superior single plants with stable yield and satisfactory disease resistance were selected as potential superior single plants. In this study, we comprehensively evaluated the economic characteristics of the mature fruit of 45 single plants based on PCA and gray correlation analysis. We ranked these plants according to their economic performance and then selected the plants with satisfactory performance as excellent germplasms. The results of this study provide a foundation for the selection of excellent cultivars and crossbreeding of C. weiningensis to achieve the ultimate goal of creating new germplasms. Additionally, the results expand the available biological data on C. weiningensis.

Experimental site and materials
All 45 superior C. weiningensis plants in this study are maintained at the C. oleifera Research Base of Guizhou University, Weining County, Bijie, Guizhou, China (long. 104˚7 0 53 0 'E, lat. 27˚11 0 53 0 'N). The average altitude of the base is 2,200 m. In this region, the daily temperature difference is large, the average annual rainfall is 909 mm, the average frost-free period is 180 d, the average annual sunshine hours is 1,812 h, and the average annual temperature is 11.5˚C. The soil at the experimental site is slightly acidic, which is suitable for the growth of C. weiningensis.
The criteria for the selection of the 45 superior C. weiningensis plants included healthy growth, a high seed-setting rate after natural pollination, a stable yield, and no diseases or insect pests. The plants that grew stably within three consecutive years were initially selected, and their height ranged from 1.5 m to 3.3 m, with a crown width of 140×160-265×270 cm. The plants were randomly designated as "WY-[serial number]". In the fruit-ripening season (August 2020), all mature fruits were collected to measure the yield per plant. Among the mature fruits, 30 were randomly selected (identified based on slight cracking of the pericarp) were bagged and labeled (Fig 1). With 24 h, measurement of the phenotypic data and fresh seed weight were completed. After measurement, the fruits were subjected to sun drying. After all samples are collected, subsequent treatment, such as drying, was performed.

Study design
After fruit collection, the fruit phenotypic and oil quality characteristics were determined and descriptive statistics and characteristic variation analysis were performed. Comprehensive assessments of the 15 economic characteristics of the mature fruit from each of the C. weiningensis plant were performed using PCA and gray correction analysis. The plants that ranked among the top ten according to both analyses were finally identified for the purpose of popularization.

Determination of economic characteristics
Six sample fruits were randomly selected from each plant. For transverse diameter, two sides of each fruit were measured, and for pericarp thickness, measurement was performed at four points of each fruit. For the remaining characteristics, measurement was performed once. The experiment was repeated six times, and the averages were obtained for each characteristic.
(1) Longitudinal and transverse diameters. The longitudinal and transverse diameters and pericarp thickness of the fruit were measured with an electronic digital caliper (Mitutoyo 500-197-20; Japan; precision, 0.01 mm). The fruit shape index was calculated as follows: Fruit shape index = transverse diameter/longitudinal diameter of the fruit.

PLOS ONE
Economic characteristics of C. weiningensis (3) Fresh seed rate. Fresh seeds were removed from the weighed fresh fruit and then weighed with an electronic balance. The fresh seed rate was calculated as follows: Fresh seed rate = fresh weight of the seeds/ fresh weight of the fruit×100% (4) Dry seed rate. Fresh seeds were dried to a constant weight and then weighed with an electronic balance, and the calculation equation was as follows: Dry seed rate = weight of the dry seeds/weight of the fresh seeds×100% (5) Dry kernel rate. Kernels were removed from the dried seeds, dried to constant weight, and then weighed with an electronic balance. The calculation equation was as follows: Dry kernel rate = weight of the dry kernels/weight of the dry seeds×100%. (6) Crown yield. The crown size of each plant was measured, and then the yield of the plant was divided by the size to obtain the crown yield (kg/m 2 ).

Oil extraction
Oil was extracted using the Soxhlet extractor method [28]. Specifically, kernels were dried at 40˚C for 24 h and then ground. Ten grams (raw weight) of each sample was wrapped in dried 12 cm filter paper (weight, W 0 ) and then dried to constant weight. After cooling, the filter paper and the sample were accurately weighed to obtain the total weight W 1 . The filter paper package with the sample was placed into a Soxhlet extractor for circular extraction for 10 h with petroleum ether. The package was then removed, dried to a constant weight and then accurately weighed. The total weight of the filter paper and residue W 2 was obtained. The oil yield was calculated as follows:

Fatty acid composition and content determination
The procedures were as follows: 1. An oil sample (100 mg, accurately weighed) was added to a conical flask.
2. Methanol (40 mL) and 2 mL of 0.5 mol/L potassium hydroxide-methanol solution were added to the oil sample. The sample underwent reflux at 75˚C. (During reflux, the flask was constantly shaken until the oil disappeared and the solution became transparent.) 3. The conical flask was removed and then cooled in chilled water to room temperature.
4. The solution in the flask was transferred into a separating funnel. The flask was washed with 20 mL of n-heptane, and the washing liquid was transferred into the funnel. Then, 20 mL of chromatography-grade pure water was added to the funnel. The solution was mixed well and then allowed to stand for stratification. (The upper layer was the ester layer, and the lower layer was the aqueous layer.) The lower layer was extracted with 20 mL of n-heptane. The extracts were combined and then washed with pure water 2-3 times until the discarded water appeared neutral. The ester layer was isolated.
5. An appropriate volume of anhydrous sodium sulfate was added to the ester solution. After filtering, the drying agent was washed with a small volume of n-heptane. The washing liquid and ester solution were merged. Afterwards, the volume of the obtained solution was increased with n-heptane to 100 mL for later use. The conditions for gas chromatography were as follows: detection, flame ionization detection (FID); column type and specification, SP2340 and 60 m × 0.25 mID 111/11 ×0.2 μm; temperature program, 50˚C (initial temperature) for 2 min followed by 170˚C (10˚C/min) for 10 min, 180˚C (2˚C/min) for 10 min and 220˚C (4˚C/min) for 22 min; inlet temperature, 250˚C; and detector temperature, 300˚C. Nitrogen was used as the carrier gas, with a separating ratio of 1:50 and a sample feed of 1 μL.
The components of fatty acids were determined based on comparisons with the retention times of samples of different fatty acid standards. The relative content of each component was calculated using the area normalization method. The experiment was repeated three times, and an average value of each component was obtained [25].

Statistical analysis
Averages, maximum values, minimum values, means, standard deviations, variation coefficients, and the matrix of comprehensive scores were calculated with Excel 2019. SPSS20.0 was used for PCA, Origin for correlation analysis and Photoshop 2019 for image processing. Gray correlation analysis was performed to comprehensively evaluate the fruit characters of different excellent plants according to the observed values of typical indices, and the gray correlation coefficient was calculated as follows: where ρ represents the resolution coefficient (set at 0.5 here). The correlation coefficient of each index of the plants was obtained. For each plant, the correlation coefficients of the characters were summed up, and an average was obtained, which was used as the gray correlation degree (r i ) of the plant. All plants were ranked according to their gray correlation degrees, and a higher r i value indicated more excellent fruit characteristics.

PLOS ONE
Economic characteristics of C. weiningensis

Variation analysis of the fruit and seed characteristics
The fruit characters of the 45 superior plants are summarized in Table 1. The genetic variation coefficients of the 15 characteristics of the plants ranged from 2.24% to 22.70%. The genetic variation coefficients of fruit diameter, fruit height, fruit shape index, fresh seed rate, dry seed rate, dry kernel rate, palmitic acid, oleic acid and arachidonic acid were small, which all fell between 0% and 15%, while single fruit weight, pericarp thickness, oil content, linoleic acid and linolenic acid had great genetic variations (all between 15% and 30%). Single fruit weight showed the greatest variation, with a range of 21.34 g and a genetic variation coefficient reaching up to 22.70%. The genetic variation coefficient of oil content was also large, which was 13.11%. These results indicated that both single fruit weight and oil content had abundant hereditary characteristics. Linoleic acid, stearic acid and pericarp thickness followed single fruit weight and oil content, whose variation coefficients were 8.40%, 17.75% and 17.36%, respectively. In contrast, the variation coefficients of oleic acid, fruit diameter and arachidonic acid were low, which were 2.24%, 6.91% and 9.10%, indicating that these characters were largely stable. From the perspective of skewness coefficient, fruit diameter, fruit length, single fruit weight, pericarp thickness, fruit shape, dry kernel yield, palmitic acid, stearic acid, linoleic acid, linolenic acid and arachidonic acid all had positive values, which indicated that their distribution curves deviated to the left (small diameter; the deviation degrees were ranked from top to bottom as follows: Single fruit weight > longitudinal diameter > stearic acid > linolenic acid > peel thickness > arachidonic acid > palmitic acid > dry kernel yield > linoleic acid > fruit diameter > fruit shape). In terms of kurtosis coefficient, the coefficients of all these characters but palmitic acid and linoleic acid were positive. This finding indicated that these characters had concentrated distribution, whose distribution curves were steeper than a normal distribution, with a sharp peak appearance. Fruit length had the greatest kurtosis coefficient, whose distribution curve were therefore the steepest. The kurtosis coefficients of palmitic acid and linoleic acid were negative, with a relatively discrete character distribution. The skewness coefficients of oil content, fresh seed yield, dry seed yield and oleic acid were also negative. These results indicated that the distribution curves of these characters were rightbiased. The kurtosis coefficients of fresh seed rate and dry seed rate were positive, whose distributions were relatively concentrated.

Correlations among characteristics
The results of the correlation analysis of the 15 characteristics were shown in Fig 2. Both diameter and height had very significant positive correlations with fruit weight. Fruit weight had a very significant negative correlation with dry kernel rate and a significant negative correlation with dry seed rate, and pericarp thickness had a very significant negative correlation, which indicated that a heavier fruit with a thicker pericarp had a lower seed rate.
Palmitic acid content exhibited very significant negative correlations with the contents of stearic acid and oleic acid. The content of this fatty acid also had a very significant positive correlation with linoleic acid content and a significant positive correlation with linolenic acid content. Stearic acid content had a very significant positive correlation with oleic acid content and very significant negative correlations with linoleic acid and linolenic acid contents. Oleic acid content showed very significant negative correlations with linoleic acid and linolenic acid contents, and linoleic acid and linolenic acid contents showed a very significant positive correlation with each other. The remarkable correlations among the components of fatty acids suggested that the content of a certain component could be evaluated based on that of another component. In addition, the degrees of correlation among the characteristics were different, which indicated information overlap among the investigated characteristics; therefore, PCA was required to reduce the number of variables.

PCA
PCA was performed on the 15 fruit character indices of the 45 C. weiningensis plants ( Table 2). Based on the criteria of cumulative contribution rate (CCR) �80% and eigenvalue (E) >1, six principal components were extracted. The CCR of the six components reached 83.06%, which meant that these components contained most of the information of the original 15 indices. These six principal components were then used as the comprehensive assessment indices for the purpose of dimension reduction.
The loading matrix of principal components can reflect the contribution of each characteristic to the principal components; moreover, it can reflect the main information conveyed by each principal component. As shown in Table 2, the first principal component contained the information about palmitic acid, stearic acid, oleic acid and linoleic acid, with a variance contribution rate of 23.586%. Oleic acid showed the most significant correlation with the first

PLOS ONE
Economic characteristics of C. weiningensis principal component, whose factor load is -0.843. The second principal component contained the information about fruit diameter, fruit shape and dry kernel rate, with a variance contribution rate of 16.584%. Fruit diameter showed the most significant correlation with the second principal component, with a factor load of -0.738. The third principal component contained the information about fruit length and single fruit weight, with a variance contribution rate of 15.551%. Fruit length showed more significant correlation with the third principal component, with a factor load of 0.719. The fourth principal component contained the information about pericarp thickness, oil content and fresh seed rate, with a variance contribution rate of 11.513%. Pericarp thickness showed the most significant correlation with the fourth principal component, with a factor load of 0.757. The fifth principal component contained the information about linolenic acid and arachidonic acid, with a variance contribution rate of 8.915%. Arachidonic acid showed a more significant correlation with the fifth principal component, with a factor load of 0.804. The sixth principal component contained the information of dry seed rate, whose factor load was 0.457. Finally, the six indices, i.e., oleic acid, fruit diameter, fruit length, pericarp thickness, arachidonic acid and dry seed rate, were finally taken as the representative indices to comprehensively assess the C. weiningensis trees.

Comprehensive evaluation based on PCA and gray correlation analysis
Data of the 45 plants were standardized, and their scores on the six principal components were obtained (Table 3). These principal components reflected the comprehensive conditions of C. weiningensis fruit at different levels, which could not be fulfilled by any single principal component. Based on the PCA scores, with the variance contribution rate corresponded by each principal component as the weight, a comprehensive evaluation function was constructed, and  the function model Z was as follows: The comprehensive scores (F) of the economic characteristics of the fruit of 45 superior C. weiningensis trees were finally obtained, according to which the trees were ranked and superior plants were screened.
In addition, based on the 15 indices of C. weiningensis, the correlation coefficients of each index were calculated. The correlation coefficients of the indices were accumulated and an average was obtained. The trees were ranked according to the averages they obtained for comprehensive evaluation.
The rankings of the comprehensive quality of the fruit of the 45 C. weiningensis plants based on PCA and gray correlation analysis are summarized in Table 4. According to the PCA, the top ten plants were ordered as follows: WY-33>WY-8>WY-1>WY-27>WY-30>WY-25>WY-38>WY-44>WY-6>WY-35. According to the gray correlation analysis, the order of the top ten plants was WY-1>WY-8>WY-38>WY-25>WY-30>WY-33>WY-6>WY-27>WY-35> WY-19. The rankings based on the two analytical methods were subjected to correlation analysis, and the correlation coefficient was 0.535, showing a very significantly positive correlation (P<0.01). Therefore, the evaluation results based on these two methods had satisfactory consistency and thus verified each other. Although the same plant might exhibit a difference in the rankings according to different evaluation methods, this difference was caused by the different weights assigned to each index during analyses. Additionally, the overall rankings of the plants based on the two methods were comparable, which indicated overall consistent evaluation trends according to the two methods. Based on these results, nine C. weiningensis plants whose fruit characters ranked in the top ten, according to both evaluation methods, were selected for promotion, which included WY-1, WY-  Table 5.
The obtained outcomes based on the PCA and gray correlation analysis were compared with the actually measured data. WY-33, which ranked first according to the PCA, ranked the first both in fruit length and single fruit weight among the plants, and WY-1, which occupied the first position according to the gray correlation analysis, ranked the first in pericarp thickness and its dry kernel rate nearly reached the highest value among all investigated plants. WN-8, which ranked the second according to both methods, ranked the first in palmitic acid, linoleic acid and arachidonic acid. The nine germplasm plants selected based on both methods exhibited excellent performance in terms of all characters. Therefore, the outcomes obtained based on PCA and gray correlation analysis were basically consistent with the actually measured data, which indicated that these nine germplasms had the potential to be applied in actual production.

Comparison with the five-year average crown yield
The obtained outcomes based on two analytical methods in this study were compared with five-year average crown yields. As shown in Table 6, WY-33, which ranked first according to the five-year crown yield, ranked first according to the PCA and the sixth according to the gray correlation analysis. WY-1, which ranked the second according to the five-year crown yield, ranked the third according to the PCA and the first according to the gray correlation analysis. There existed inconsistencies. For instance, WY-16 ranked the 35th in terms of the crown yield and the 44th according to both analytical methods. Although the overall fruit

PLOS ONE
Economic characteristics of C. weiningensis quality of this plant was relatively inferior, its performance in yield was not so poor. Also, although some plants had comparatively low yields, their overall fruit quality was relatively satisfactory (e.g., WY-28 and WY-29). Despite the inconsistencies, the assessment outcomes based on PCA and gray correlation analysis were consistent with the actual yields to a certain extent, which indicated the reliability of the assessment outcomes in this study. In addition, based on the results of the comparison with the five-year average crown yield, the PCA outperformed the gray correlation analysis.

Discussion
The fatty acids of edible vegetable oil are mainly composed of monounsaturated fatty acids (MUFAs) and polyunsaturated fatty acids (PUFAs

PLOS ONE
Economic characteristics of C. weiningensis Notes: Data on crown yield are presented as the mean ± standard deviation. PCA, principal component analysis; GCA, gray correlation analysis. than saturated fatty acids benefits blood sugar and blood lipid control, which can partially reduce the risks of hyperinsulinemia, pancreatitis and β-cell dysfunction induced by excessive intake of saturated fatty acids and the risk of invasion of renal cell cancer to renal cells [30]. Presently, cardiovascular and cerebrovascular diseases and diabetes have high incidence; therefore, C. weiningensis oil has great potential for promotion. Furthermore, with the improvement of people's living standards, the demand for edible vegetable oil is increasing. China has become a large consumer and importing country of edible vegetable oil; approximately 70% of the oil consumed in the country is imported. At present, edible vegetable oil in China is mainly extracted from rapeseed, peanut, soybean and other herbal oil crops. The country is unlikely to develop woody oil production by sacrificing cultivated land. However, developing woody oil crops in mountainous areas may be a feasible solution to this issue. This solution may not only reduce the current reliance of China on import for edible oil but also conserve areas of cultivated land for grain crops, thereby maintaining food security in China [31]. Current methods adopted for the comprehensive assessment of fruit economic characteristics mainly include PCA, gray correlation analysis, DTOPSIS and the rationalization-satisfaction method. PCA and grey correlation degree are both commonly used multivariate statistical methods. Although their analysis principles are different, they can both be used to comprehensively and objectively evaluate the multiple characters of fruit. PCA offers advantages in that it does not entail artificial weighting of each characteristic and the contribution rates of the principal components function as the weights used in other evaluation methods; these features enable the entire assessment process to be objective. Gray correlation analysis can solve the difficulty in index selection when the indices are in an excessively large number, and thus realizes comprehensive, all-round evaluation. Based on the final comprehensive scores and ordering outcomes of the superior plants in this study, the superior C. weiningensis plants tended to be those with heavy single fruits, a high kernel rate and high palmitic acid, oleic acid, linoleic acid and arachidonic acid contents.
Oil yield and fatty acid composition are important indices for the comprehensive assessment of the economic characteristics of C. oleifera. The oil yield and fatty acid composition of C. weiningensis differ to various degrees from those of Hainan C. oleifera [32]: Both the average oil yield and average oleic acid content of C. weiningensis are higher than those of C. oleifera. Moreover, the variation coefficients of their characteristics also differ. The variation coefficients of C. weiningensis are generally lower than 0.1%, whereas in Hainan C. oleifera, those of linoleic acid (46.12%), oleic acid (6.39%), palmitic acid (13.41%), stearic acid (21.53%) and linolenic acid contents (39.7%) are remarkable. Growth regions and cultivar types have tremendous impacts on the properties of tea seed oil, such as oil components and fatty acid composition [33]. Such impacts are mainly caused by precipitation, illumination, temperature, altitude and topography [34][35][36][37][38]. Weining County has a high average altitude, and the oleic acid, linoleic acid and linolenic acid contents of C. weiningensis were 78. 35 [42] (with all the above-listed regions having average altitudes far lower than the Weining region). Camellia weiningensis has a higher oleic acid content than C. oleifera in 80% of these regions, a lower linoleic acid content than C. oleifera in 75% of these regions and the second highest linolenic acid (following that of C. oleifera in Liuyang, Hunan). The comparative results did not show the particular advantage of C. weiningensis due to its high cultivation altitude. The content of saturated fatty acids in tea oil tends to increase with increasing latitude of the planting area, whereas that of unsaturated fatty acids exhibits the opposite trend [43,44]. Between the altitudes of 297 m and 880 m, the oil content and chemical components of tea oil of wild C. oleifera vary greatly; however, the variations do not have noticeable associations with altitude [45]. Between altitudes of 1,900 and 2,500 m, the unsaturated fatty acid content of walnut oil shows a wave-like fluctuations but altitude has no noticeable influence on the oil's linolenic acid content [46]. In green coffee beans, the four principal fatty acids, i.e., palmitic acid, linolenic acid, stearic acid and oleic acid, show moderately negative correlations with altitude [47]. Therefore, although the accumulation of oleic, linoleic and linolenic acids during seed development of C. oleifera varies greatly among altitudes, the differences are mainly attributed to differences in climatic factors, and the fatty acid composition of seeds is not significantly correlated with altitude. Thus, altitude cannot be used as an evaluation criterion for the fatty acid components of C. oleifera seeds.
The mechanisms of action of site conditions on the oil yield and fatty acid composition of C. oleifera have been reported. Soil organic matter (SOM) may serve as a limiting factor for oil yield, whereas available phosphorus (AP) deficiency is the main limiting factor for oil yield in the central subtropical regions of China [38]. During seed development, the relative activity of fatty acid desaturase 2 (FAD2), which is associated with the syntheses of oleic acid and linoleic acid, continuously decreases, whereas that of FAD3 gradually increases. The differences in the relative activity of FAD2 among regions of different altitude might be associated with variations in the average daily temperature and average daily temperature difference; suitable temperature, a large temperature difference and low precipitation all promote oil accumulation in seeds, thereby increasing the oil content [48,49]. According to the literature, the net photosynthetic rate has a very significant correlation with oil content (P<0.001) and a significant correlation with the kernel rate (P<0.01). Acid rain significantly enhances the leaf photosynthetic rate of woody tree species [50]. In Weining, acid rain occurs frequently, which might constitute one cause of the higher oil yield of C. weiningensis than of C. oleifera in half the regions mentioned above. However, in this study, we did not explore the mechanism underlying the high oil yield of C. weiningensis, which may be an important direction for future research. Analysis of the mechanism of action each factor will undoubtedly benefit the improvement of site conditions for C. weiningensis, thereby enhancing the yield and quality of the tea oil.
In tea seed oil, direct transformation processes occur from palmitic acid to stearic acid, from stearic acid to oleic acid, from oleic acid to linoleic acid and from linoleic acid to linolenic acid, and the selection of a suitable fruit harvesting time can improve the proportion of unsaturated fatty acids in tea seed oil [48]). Therefore, the selection of the optimal fruit cropping time can bring the oil, oleic acid, linoleic acid, and linolenic acid contents and fresh seed and kernel rates under balanced control [51]. Late August to early September is a suitable harvest time for C. oleifera in Honghua, Zhejiang [52]. With the maturation of C. oleifera seeds, the unsaturated fatty acid content in seed oil increases and the levels of antioxidants, such as vitamin E and βsitosterol, increase; therefore, the time after October and before fruit fall is the best fruit-harvesting time for C. oleifera [53].
In this study, we comprehensively assessed 45 superior C. weiningensis cultivars using PCA and gray correlation analysis and screened nine superior cultivars for breeding. In addition, the coefficients of genetic variation for oleic acid, transverse diameter and arachidonic acid were all lower than 10%, which has not been reported before. Presumably, the low levels of variation in these components are due to SOM content, harvest time, high altitude and frequent acid rain. However, to verify this presumption, additional systematic research needs to be conducted in the future. In addition, currently, the popularization of C. weiningensis is at the superior variety selection stage, and mass production has not yet begun. Therefore, other indices to evaluate the superiority of the selected varieties, such as unit output, cannot yet be calculated; this topic constitutes another important direction for further research.